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ABSTRACT 

In this paper we use grids defined in Cartesian coordinates in place of grids defined in polar and spherical polar 
coordinates to solve an elliptic partial differential equation in two and three space dimensions. This transforming of the 
grids makes using both the interpolation and restriction operators more simple and moreover they give convergence rate 

better than operators defined in polar and spherical polar coordinates. A full multigrid method with W ( y 1 , Yi ) cycle 

with line relaxation and both restriction and interpolation operators defined in Cartesian form are used. Finally numerical 
examples have been given. 

KEYWORDS: Multigrid Method, Numerical Analysis, Elliptic PDE, Cartesian, Polar and Spherical Polar Coordinate 
Systems 

1 INTRODUCTION 

Consider the two dimensional Poisson's equation defined in polar coordinate system: 

U rr +-U,+± r U eg =-F(r ,0)en 
r r 

subject to U = G on d£2 where 3/2 is a boundary of a quarter of a unit circle defined by: 

n={(r,0): 0 <r <1, 0 < 0 <—} (1.1) 

2 

and Consider the two dimensional general linear elliptic partial differential equation defined in polar coordinate 

system: 

(1 + sin 0) U rr + — (1 + sindcosd) U r + — - sin dcos 9) U eg 

1 ' ■ .2 n 2 n\ tj , / 2 



+ —(sin 2 8 - cos' 6) U e +-(cos z 6 - sin 2 6) U r0 = - F(r,8) e Q 

subject to U — G on 3i2 where 3/2 is a boundary of a quarter of a unit circle defined by 



n={(r,0): 0 <r <1,O<0<-} (1.2) 

2 

also consider three dimensional Poisson's equation in spherical polar coordinate system: 
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U rr +-U r +-U ff +--^-U f + 2 \ 2 U ee =-F(r,<p,0)en 
r r 2 wv r 2 sin <p w r sin <p 

subject to U = G on 9i2 

where dn is a boundary of a eighth of a unit ball defined by (1.3) 

n={(r, q>, 6) : 0 < r <1, 0 <(p, 6<—} 

2 MULTIGRID METHOD AND TRANSFORMATION OF THE GRIDS 

For the multigrid method, we needs a sequence of grids, then replacing each term in equation 1.1 by the 
corresponding finite-difference approximation in the r 6 -plane see [1] and [2]. In this paper we will solve these problems 
by using transformed grids defined in Cartesian coordinate system to overcome the difficulties in both interpolation and 
restriction operators. In the two dimensional space the difficulties of the interpolation operator is how to determine the 
interpolating values for points lying on the curved boundary and also for points lying in the neighborhood of the center of 
the circle of the fine grid and the difficulties of the restriction operator is how to calculate weights of the full-weight 
restriction operator see [1]. In the three dimensional space (in addition to the above difficulties of the interpolation 
operator) the difficulties of the interpolation operator is how to determine the interpolating values for points lying on the 
curved surfaces and also for points lying in the neighborhood of the lines passing through the center of the sphere and any 
of the poles of the sphere. In the transformed grids defined in Cartesian coordinate system we will use the a full multigrid 

method with W ( YnYi) cycle, full weight residual restriction operator is used to get the residuals at coarse-grid points 

from the residuals computed at points on the fine grids. Also interpolation operator is used to get the approximation 
solution at fine-grid points from the solution computed at the coarse-grid points and even-odd line Gauss-Seidel relaxation 
is used. Numerical examples have been given. 

3 FINITE DIFFERENCE DISCRETIZATON IN POLAR FORM AND IN SPHERICAL POLAR FORM 

We will use the following finite difference approximation in the two dimensional space in polar form: 

U =h~ 2 ( u . , . +u. , . —2u. ■ ) 

U ee =(d0y 2 (u i . +l+ u iJ _ l -2u i .) (3i) 

U e =(2S0)- 1 (u iJ+1 -u iJ _ 1 ) ^ 

U rg =(4hddy l (u i+l j+1 +u i _ 1J _ 1 -u i+1J _j -u,^ J+1 ) 
F(r i ,9 i )=f iJ 



Impact Factor (JCC): 4.2949 
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form: 



r - ih, i = 1, 2, N 

and we will use the following finite difference approximation in the three dimensional space in spherical polar 



U rr = h~ 2 (u i+ljk +u i _ 1Jik -2u i j k ) 

U r =(2h)- 1 (u i+1Jik -u i _ 1JA ) 

U„ =( S<p)- 2 (u iJ+l k +u iJ _ l k -2u iJ k ) 

U (p =(2 8(p)- 1 (u i j+l k ~u-j_ l k ) 

U ge =( 80y 2 (u i j k+1 +u iJ k _ 1 -2u iJ k ) 

F(r i , <Pj ,e k )=f i>j>k 

r = ih, i = 1, 2, N 



(3.2) 



Substitute each term of equation 3.1 in equation 1.1 we get the following system of equations in a matrix form 



given by: 



-A 


B 


0 




0 


B 


-A 


B 




0 


0 


B 


-A 


B 


0 










0 


0 




B 


-A 


B 


0 




0 


B 


-A 









u 2 






}*N-1 _ 




b N-I 



(3.3) 



where 1 <i , j < N -1 , h = — , dd = 

N +1 2(N+1) 



are the mesh size and 



A =2 



1 +• 



(isey 



(3.4) 



B 



(isey 



where 1 <i , j < N —1 , U -j is an approximation to the exact solution U ( r i , 0- ), b - is the right hand 



/ 1 

side, h = , 60 = 

N +1 2(N+1) 



, are the mesh sizes, where is the number of interior points in each 



direction. Substitute each term of equation 3.1 in equation 1.2 we get the following system of equations in a 
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matrix form given by: 



-A 




0 




0 


B 2 


-A 


B, 




0 


0 


B 2 


-A 


B> 


0 










0 


0 




B 2 


-A 


B> 


0 




0 


B 2 


-A 



u 1 




~*1 ' 


u 2 










Pn-1. 



(3.5) 



where 1< i, j <N -1, h = / r , 60 = 1 



N +1 2(N+1) 



are the mesh size and 



A = 2 



B,= 



1 + —sin2 0 + — - — T (l -—sin 2 0) 
2 (iS0) 2 2 



1 1 . cos 2 0 
T (l sin 2 0) — 



1 i . „„. cos 2 0 

— — -^(7 sin 2 (9) + — 7— 



(3.6) 



where 1 <i , j <N —1 , U •• is an approximation to the exact solution U( r., 0. ), b •• is the right hand 

, are the mesh sizes, where JV is the number of interior points 



side, = , 60 — 

N +1 2(N+1) 



Also substitute each term of equation 3.2 in equation 1.3 we get the following system of equations in a matrix form given 
by: 



-A 


B 


0 




0 


B 


-A 


B 




0 


0 


B 


-A 


B 


0 










0 


0 




B 


-A 


B 


0 




0 


B 


-A 



u 1 






u 2 




b 2 






Ps-1_ 



where,/ < i,j ,k <N -1, U iJ>k 
U(r iy (p., 0 k ), b ijik is the right hand side, ,tl 



A =2 



1 +■ 



(iSq>) 2 (i Sdsin <p) 2 
Impact Factor (JCC): 4.2949 



(3.7) 

is an approximation to the exact solution 

, 6(0 = 60 = are the mesh sizes, and 

N +1 * 2(N+1) 

(3.8) 
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B = 



( i 56 sin <p) 2 

where N is the number of interior points in each direction. 

4 TRANSFORMATION OF THE GRIDS FROM A POLAR AND A SPHERICAL POLAR COORDINATES 
TO A CARTESIAN COORDINATE SYSTEM 

The two dimensional domain defined in polar coordinates given by 
n={(r,0): 0 <r <1, 0<0<-} 
may be transformed into a domain defined in a Cartesian coordinates given by: 

n={(x,y ) : 0 <x,y <lj 



2h 




r=0 2h 0 = o 



TIT 



0 = 0 












7 


a 


<> 




4 


5 


6 




i 


2 


3 





Figure 1: 2D-Transformation from a Polar to a Cartesian Grid 

(see Figure 1) 

also the three dimensional domain defined in spherical polar coordinates given by 
n={r ,(p,6) : 0 <r <1, 0 <<p, 0 <—} 

may be transformed into a domain defined in a Cartesian coordinates given by: 

n={(x,y,z): 0 <x,y,z <1} 
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Figure 2: 3D-Transformation from a Polar to a Cartesian Grid 



(see Figure 2) 



5 FULL MULTIGRID METHOD 



The elliptic equation have been solved using full multigrid method with W ( Yi > Yi ) -cycle that has a sequence 



1 



of 5 grids f2 kh have been used with mesh sizes kh, k=2",«=0,l,...4, where h = ^~> w i tn tne coarse grid correction 
(CGC) and has the following components: see [2], [3] and see Figure 3 



1. W( Yj, Yi) Cycle with: 

a) Vj is the number of relaxation sweeps before coarse grid correction (CGC). 

b) Coarse grid correction (CGC). 

c) V 2 is the number of relaxation sweeps after coarse grid correction (CGC). 

2. Red-Black Line Gauss-Seidel Relaxation: The most efficient smoothing iteration (relaxation) process is the red-black 
Gauss-Seidel iteration for lines which gives better results for polar problems than the point relaxation: 

• For the two dimensional space the general grid can be defined by: 




Figure 3: FMG W ( YnYi) c y cle 



Q k = {(vh k , juh k ), v, jueZ, Z is the set of integer numbers 



(5.1) 



then split Q k into red (even) lines: 




e Q k , v is even 



(5.2) 
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and the remaining lines are black (odd): 

Q k ={(vh k , juh k ), e& k , vis odd] (5.3) 

• For the three dimensional space the general grid can be defined by: 

£2 k ={( vh k , juh k , yh k ), v, ju, ysZ, Z is the set of integer numbers} (5.4) 

then split £2 k into red (even) lines: 

A" ={ V K > J^k > fan ) G Q k ,: ( V + JU) is even} (5.5) 
and the remaining lines are black (odd): 

Q\ ={( vh k ,juh k ,yh k )e n k ,:(v + ju) is odd} (5.6) 

First the solutions are calculated at the points of the even (odd) lines using the solution at the points of the odd 
(even) lines and vise-versa, so for each line we have to use Gauss elimination method to solve a system of linear equation 
for equation 3.5 and for equation 3.7. 

t h 2 h 

3. For Coarse to Fine Interpolation with Operator 1 : It takes coarse-grid vector (solution or correction) U defined on 
the coarse grid f2 2 h and produces fine grid vector U h defined on the fine grid £2 h according to the rule 
T 2h l 



1 2 h U 2h —U h . the two and three dimensional space, we have the following components of U h 



* For the two dimensional space the components of U h are as follows: 
h 2h 

U 2 i2j j for common points in the two grids (point 5) (5.7a) 

h 1 / 2h 2h \ 

**2i+l 2i ~ — ' / i) f° r points in the r = x direction, (points 4, 6) (5.7b) 

^2i2j+l ~ +7 ) for points in the 0—y direction, (points 2, 8) (5.7c) 

h 1 / 2h 2h 2h 2h \ 

for intermediate points, (points 1, 3, 5, 7) 
see Figure 1 

h 

* For the three dimensional space the components of U are as follows 

U 2i 2j 2k =u i'j k ^ or common points in the two grids (5.8a) 
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u 2i+l,2j,2k = ~^( u iJ,k +u i+l,j,k ) for points in the r = x direction (5.8b) 

u 2i,2j+l,2k =y( M , 2 j,/t +u i,j+l,k ) for points in the (p=y direction (5.8c) 

u 2i,2j,2k+l = Y^ U fj,k +u i H j,k+l ) for points in the ^ = Z direction (5.8d) 

U 2i+l,2j+l,2k = j( U iJ,k +U i+l,j,k +U iJ+l,k +U i+l,j+l,k ) (5 - 8 ^ 
for points in the r (p — Xy plane direction 

u 2i+l,2j,2k+l = ^( u iJ,k +u i+l,j,k +u iJ,k+l +u i+l,j,k+l) ( 58f ) 



for points in the V 0 — XZ plane direction 

«ti,2J + l,2k + l =j( ^Ik +<W + U ? h j,k + 1 + U i,j + l,k + l ) ^ 

for points in the (p0 — yz plane direction 

h 1 / 2h 2h 2h 2h 

U 2i+l,2j+l,2k+l ~~^( U i,j,k +U i+l,j,k +U i J +1 ,k +U i J ,k +1 

+ u i+l,j+i,k +u i+i,j,k+i +u i,j+l,k+i +u i +1 ,j +1 ,k +1 ) (5.8h) 

For intermediate points see Figure 2 

4 Fine to Coarse Restriction Operator: The restriction operator denoted by I £ h takes the residual vector R h computed on 

the fine-grid and transfer it to the coarse-grid according to the rule: I^ h R h = R 2h , where the components of R(j of 

R 2h are given by two restriction operators: half-weight and full-weight restriction operator; since we use the line 
relaxation, then the full weight operator will be more effective than the half- weight operator see [2] and [3]: 

• For the two dimensional space the full weight restriction operator R 2h will be defined by: 

R i,j =J^[ R 2i-l,2j-l + R 2i-l,2j+l + R 2i+l,2j-l +R 2i+l,2j+l + ( 5 - 9 ) 
2 ( R 2i,2j-l +R 2i,2j + 1 +R 2i-l,2j +R 2i + l,2j ) +4R 2i,2j] 
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The matrix form of the full-weight operator 1 h takes the form 



j 2h _ ± 

* ~16 



12 1 
2 4 2 
12 1 



2h 



For the three dimensional space the full weight restriction operator R ' will be defined by: 
R i,j,k = Y R 2i,2j,2k + J S 1 + V S 2 + 1T S 3 



where 



C _ n h in* 4- J? fe 4- J? A 

J 1 ~ "2i-l,2j ,2k " r "2i+l,lj ,2k ~*~ " 2 i ,2 j -1 ,2 k " r A 2i ,2j +1 ,2k 

, 7) A . nh 

~ r "2i,2j ,2k-l ~ r "2i,2j+l,2k+l 

p d A 4- J? A 4- J? * 4- J? * 

° 2 ~ "2i-l,2j-l,2k "2i+l,2j -1,2k " 2 i -1 ,2 j +1 ,2k ~*~ " 2 i +1 ,2 j +1 ,2 k 

i D A 4-Ti h 4-J?' 1 4-J?' 1 

' r "2i,2j-l,2k-l ' r "2i,2j+l,2k-l ~ " 2 i ,2 j -1 ,2k +1 ~ rly 2i,2j- 

ijfh _|_ n /; _|_ r> /i , n /i 

^^21-1,2] ,2k -1 2 i +1 ,2j ,2 k -1 ^ 2 i -1 ,2j ,2 k +1 ~*~ 2i +1 ,2j ,2k +1 

P _pS 4-J?' 1 4- J? * 4-J? h 

J i — "-21-1,2 j -1,2k ~*~"2i+l,2j-l,2k-l ~*~ " 2 i -1 ,2 j +1 ,2 k -1 ~*~ " 2 i -1 ,2 j -1 ,2 k +1 

_i_ JO A i n' 4-J?' 1 4-J?' 1 

' TJV 2i-l,2j+l,2k+l " l " A 2f+I,2i-I,2il+i " 1 " ^2 i +7,2./ +7,2/fc -7 2 / +7,2 j +7 ,2 A: +7 



;+7,2t+7 



The matrix form of the full-weight operator takes the form: 



j 2h _ 



/ 


2 


/" 




"2 


4 


2" 




2 




2 


+ 




S 




+ 


/ 


2 


i 




2 


4 


2 





i5 



6 TEST PROBLEMS, RESULTS AND PERFORMANCE 



12 1 
2 4 2 
12 1 



,2A 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



In this section we report the results and performance of some test problems All our implementations in Fortran 
were executed on a Pentium 4 PC using FORTRAN90 workstation compiler. The accuracy of the method is measured by 
both L 2 norm of both defect and error and the max norm of both defect and the error. 

6.1 Two Dimensional Poisson's Equation Defined in Polar Coordinate System 

For the two dimensional Poisson's equation defined in polar coordinate system (equation 1.1), consider the 
following test problems: 
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1 . Test Problem 1 : 

2. Test Problem 2: 

3. Test Problem 3: 



t/ , =e 



r 2 sin Gcos 0 



V exact =sin[r( 3 cos 9 + sin 9] 



U pr . =sin[ r 2 sin 9 + sin 9] 



Table 1: Fortran Implementation of Test Problem 1 Using FMW (2, 1) on PC 
Test Problem 1: Equation 6.1 



(6.1) 
(6.2) 
(6.3) 





i 

Y 
1 
4 


0.232420D- 


-07 


0.232420D- 


-07 


0.389638D- 


-01 


0.389638D- 


-01 




0.285286D- 


-05 


0.106780D- 


-05 


0.104190D- 


-01 


0.613996D- 


-02 


Equation 6.1 


1 

8 
1 

16 


0.111818D- 
0.157221D- 


-04 
-04 


0.395417D- 
0.327736D- 


-05 
-05 


0.265074D- 
0.672460D- 


-02 
-03 


0.132351D- 
0.307192D- 


-02 
-03 




1 

32 


0.9229 16D- 


-04 


0.492099D- 


-05 


0.167608D- 


-03 


0.487335D- 


-04 




1 

64 


0.303821D- 


-03 


0.146161D- 


-04 


0.408888D- 


-04 


0.183232D- 


-04 



Table 2: Fortran Implementation of Test Problem 2 Using FMW (2, 1) on PC 



Test Problem 2: Equation 6.2 







Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 






0.331 118D-05 


0.331 118D-07 


0.281695D-01 


0.281695D-01 






0.488095D-05 


0.23487 1D-05 


0.132323D-01 


0.814887D-02 


Equation 6.2 




0.112966D-03 


0.238629D-04 


0.36506 1D-02 


0.187780D-02 












16 


0.185343D-03 


0.195310D-04 


0.9322 17D-03 


0.439905D-03 




32 


0.113296D-03 


0.714172D-05 


0.234723D-03 


0.106309D-04 




1 

64 


0.501861D-04 


0.203366D-05 


0.580549D-04 


0.256458D-04 



Table 3: Fortran Implementation of Test Problem 3 Using FMW (2, 1) on PC 



Test Problem 3: Equation 6.3 



Test problem 


h 


Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 






0.479 150D-08 


0.479 150D-08 


0.139922D-01 


0.139922D-01 






0.140393D-05 


0.538168D-06 


0.356 133D-02 


0.249522D-02 


Equation 6.3 




0.875335D-05 


0.34963 1D-05 


0.892982D-03 


0.533545D-03 












16 


0.488703D-05 


0.146127D-05 


0.223801D-03 


0.123341D-03 




32 


0.302690D-05 


0.661980D-06 


0.55775 1D-04 


0.297402D-04 




1 

64 


0.149890D-05 


0.256758D-06 


0.135154D-04 


0.719024D-05 



6.2 Two Dimensional General Elliptic Partial Differential Equation Defined in Polar Coordinate System 

For the two dimensional general elliptic partial differential equation defined in polar coordinate system (equation 
1.2), consider the following test problems: 
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1 . Test Problem 4: 



2. Test Problem 5: 



exact 



U_ t =e°- 5r2sindcose +2r 2 sin Ocos 6 

r 4 (sin Ocos 0) 2 



^ exact ^ 



3. Test Problem 6: V ' exacl =r 6 (sin Ocos 0) 3 +3r 4 (sin Ocos Q) 2 +6r 2 (sin Ocos 0) 

4. Test Problem 7: U exact = s * n ( r 2 Sln ® cos & + 1 ) 

Table 4: Fortran Implementation of Test Problem 4 Using FMW (2, 1) on PC 







Test Problem 4: Equation 6.4 






Test problem 


h 


Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 






0.487803D-06 


0.487803D-06 


0.823281D-01 


0.823281D-01 






0.293 176D-04 


0.179425D-04 


0.130832D-01 


0.681229D-02 


Equation 6.4 




0.109685D-02 


0.35248 1D-03 


0.324392D-02 


0.128157D-02 














16 


0.117734D-02 


0.234885D-03 


0.942707D-03 


0.273020D-03 




32 


0.624465D-03 


0.106234D-03 


0.383377D-03 


0.806606D-04 




1 

64 


0.561402D-03 


0.452083D-04 


0.144482D-03 


0.258706D-04 



Table 5: Fortran Implementation of Test Problem 5 Using FMW (2, 1) on PC 







Test Problem 5: Equation 6.5 






Test problem 


h 


Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 




I 

2 

l 

4 


0.3401 16D-06 


0.3401 16D-06 


0.355722D-01 


0.355722D-01 




0.135265D-04 


0.700796D-05 


0.939620D-02 


0.436022D-02 


Equation 6.5 


1 

8 
1 

16 


0.738395D-04 
0.130248D-03 


0.267780D-04 
0.230902D-04 


0.237930D-02 
0.603557D-03 


0.977388D-03 
0.226562D-03 




1 

32 


0.780370D-04 


0.114822D-04 


0.159740D-03 


0.543615D-04 




1 

64 


0.175214D-03 


0.912809D-05 


0.410080D-04 


0.136098D-04 



Table 6: Fortran Implementation of Test Problem 6 Using FMW (2, 1) on PC 







Test Problem 6: Equation 6.6 






Test problem 


h 


Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 






0.814372D-07 


0.814372D-07 


0.163969D+00 


0.163969D+00 






0.418275D-04 


0.190097D-04 


0.299146D-01 


0.153482D-01 


Equation 6.6 




0.118158D-02 


0.385504D-03 


0.753403D-02 


0.318098D-02 














16 


0.138994D-02 


0.263022D-03 


0.199151D-02 


0.694873D-03 




32 


0.749245D-03 


0.119769D-03 


0.641346D-03 


0.172192D-03 




1 

64 


0.478077D-03 


0.4490 16D-04 


0.208378D-03 


0.4623 17D-04 
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Table 7: Fortran Implementation of Test Problem 7 Using FMW (2, 1) on PC 







Test Problem 7: Equation 6.7 










Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 




Equation 6.7 


i 

2 
1 

4 
1 
8 
1 

16 
1 

32 
1 

64 


0.150675D-06 
0.619946D-05 
0.457767D-04 
0.371336D-04 
0.492353D-04 
0.180558D-03 


0.150675D-06 
0.351387D-06 
0.147752D-04 
0.101259D-04 
0.546994D-05 
0.756322D-05 


0.673306D-02 
0. 14333 1D-02 
0.354350D-03 
0.926256D-04 
0.267029D-04 
0.762940D-05 


0.673306D-02 
0.783644D-03 
0.187174D-03 
0.4433 16D-04 
0.110106D-04 
0.26646 1D-05 



6.3 Three dimensional Poisson's Equation Defined in Spherical Polar Coordinate System: For the three dimensional 
Poisson's equation defined in spherical polar coordinate system: (equation 1.2), consider the following test problems: 

1. Test Problem 8: U exact = r 6 sin 4 fCOS 2 fsin 2 0COS 2 0+1 (6.8) 



2. Test Problem 9: U exact -sin( r 6 sin 4 (p cos 2 (p sin 2 Qcos 2 9 + 1) (6.9) 



3. Test Problem 10: U exact = e r 6 sin 4 <pcos 2 <psin 2 Ocos 2 0+1 (6. 10) 



Table 8: Fortran Implementation of Test Problem 8 Using FMW (2, 1) on PC 







Test Problem 8: Equation 6.8 








h 


Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 


Equation 6.8 


i 

2 
1 

4 
1 
8 
1 

16 


0.813771D-06 
0.352586D-03 
0.143917D-03 
0.144340D-02 


0.813771D-06 
0.150108D-04 
0.932530D-04 
0.230529D-03 


0.400865D-02 
0.212562D-02 
0.674367D-03 
0.133991D-03 


0.400865D-02 
0.142138D-02 
0.456735D-03 
0.142577D-03 



Table 9: Fortran Implementation of Test Problem 9 Using FMW (2, 1) on PC 







Test Problem 9: Equation 6.9 








h 






Equation 6.9 


i 

2 
1 

4~ 
1 
8 
1 

16 


0.144218D-06 
0.185749D-03 
0.796397D-04 
0.144420D-02 


0.144218D-06 
0.790830D-04 
0.518679D-04 
0.205533D-03 


0.20964 1D-02 
0.108546D-02 
0.33301 1D-03 
0.701547D-04 


0.20964 1D-02 
0.728700D-03 
0.228682D-03 
0.725672D-04 
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Table 10: Fortran Implementation of Test Problem 10 Using FMW (2, 1) on PC 







Test Problem 10: Equation 6.10 






Test problem 


h 


Max N(R) 


L 2 (R) 


Max N(E) 


L 2 (E) 


Equation 6.10 


i 

2 
1 

4 
1 

8~ 
1 

16 


0.120721D-05 
0.975 143D-03 
0.39690 1D-03 
0.489924D-02 


0.120721D-05 
0.415379D-03 
0.257905D-03 
0.709457D-03 


0.111215D-01 
0.598264D-02 
0.193477D-02 
0.377893D-03 


0.111215D-01 
0.399 128D-02 
0.130048D-02 
0.402946D-03 



7 CONCLUSIONS 

In the previous paper [1] there were some difficulties in determining both interpolation and restriction operators, 
in the two dimensional space the difficulties of the interpolation operator is how to determine the interpolating values for 
points lying on the curved boundary points and also for points lying in the neighborhood of the center of the circle of the 
fine grid and the difficulties of the restriction operator is how to calculate weights of the full-weight restriction operator see 
[1]. In the three dimensional space (in addition to the above difficulties of the interpolation operator) the difficulties of the 
interpolation operator is how to determine the interpolating values for points lying on the curved surfaces and also for 
points lying in the neighborhood of the lines passing through the center of the sphere and any of the poles of the sphere. In 
this paper we have used transformed grids defined in Cartesian coordinate system to avoid these difficulties, in addition we 
have found that the convergence rate is better than that with grids defined in polar and spherical polar coordinate system, 

this is because the round-off error depends on both (i , and (i , , 8(p, <5#) while in Cartesian coordinate system the 

round-off error is independent of these factors. 
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